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ABSTRACT 

We have investigated the oxygen and nitrogen chemical abundances in extremely compact star- 
forming galaxies (SFGs) with redshifts between ~0.11 and 0.35, popularly referred to as "green 
peas". Direct and strong- line methods sensitive to the N/O ratio applied to their Sloan Digital 
Sky Survey (SDSS) spectra reveal that these systems are genuine metal-poor galaxies, with mean 
oxygen abundances ~ 20% solar. At a given metallicity these galaxies display systematically large 
N/O ratios compared to normal galaxies, which can explain the strong difference between our 
metallicitics measurements and previous ones. While their N/O ratios follow the relation with 
stellar mass of local SFGs in the SDSS, we find that the mass-metallicity relation of the "green 
peas" is offset >0.3 dex to lower metallicities. We argue that recent interaction-induced inflow 
of gas, possibly coupled with a selective metal-rich gas loss, driven by supernova winds, may 
explain our findings and the known galaxy properties, namely high specific star formation rates, 
extreme compactness, and disturbed optical morphologies. The "green pea" galaxy properties 
seem to be not common in the nearby universe, suggesting a short and extreme stage of their 
evolution. Therefore, these galaxies may allow us to study in great detail many processes, such 
as starburst activity and chemical enrichment, under physical conditions approaching those in 
galaxies at higher redshifts. 

Subject headings: galaxies: dwarf — galaxies: abundances — galaxies: starburst — galaxies: evolution 



1. INTRODUCTION 

The cosmological relevance of local starbursts 
is due to their similarities with their analogs at 
high redshift. They may therefore provide ex- 
tremely valuable laboratories to study with high 
resolution and sensitivity relevant physical pro- 
cesses associated with the starburst activity in 
much better detail. One of the most suitable en- 



vironments for these purposes can be the intrigu- 
ing class of very compact, extremely star- forming 
galaxy (SFG) at low redshift (0.11 < z < 0.35) 
recently disco vered by volunteers in the "Galaxy 
Zoo" project (jLintott et al. I l2008l) . These galax- 
ies, popularly referred to as "green peas" (here- 
after GPs ), were first report e d and studied at some 
length bv lCardamone et al. | (I2009L hereafter C09), 
who classified more than one hundred of them. 
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The GPs are rare objects, located in lower-density 
environments and spectroscopically characterized 
by very faint continuum emission and strong op- 
tical emission lines, namely [Oiii]A5007, with an 
unusually large equivalent width of up to ^lOOOA. 
These properties along with their extreme com- 
pactness (unresolved in Sloan Digital Sky Survey 
(SDSS) imaging) are the reason for their green col- 
ors and point-like appearance in the SDSS plates. 
According to C09, GPs are low-mass galaxies (stel- 
lar masses M*~10 8 ' 5 to 10 10 M Q ) with extremely 
high star formation rates (SFRs up to 30 Mq 
yr _1 ), as derived from their emission lines and UV 
luminosity. Therefore, their specific SFRs (SSFRs, 
i.e., SFR per unit of mass) are among the highest 
inferred in the nearby universe. The GP proper- 
ties are the n consistent with th e picture of galaxy 
downsizing ( Cowie et al. 19961 ). where the sites of 



active star formation shift from high-mass galaxies 
at early times to lower mass systems at later epoch 
(see alsolBundv et al.l l2005 ICowie fc Bargeif2 008; 



Noeske et alj|2007l ). 

The motivation for this Letter was one of 
the most striking results obtained by C09: the 
nearly constant solar metallicity of the GPs 
(log[0/H]+12~8.7 in median) over more than two 
orders of magnitudes in stellar mass. These results 
would imply strong constrains to the nature of 
GPs and their evolutionary status. In this Letter 
we examine in more detail the GP chemical abun- 
dances, namely oxygen and nitrogen, and their re- 
lation with the stellar mass, using the same SDSS 
spectra but a different methodology than in C09. 
This analysis indicates instead, that GPs can be 
considered as a genuine population of metal-poor 
galaxies. The results are compared with a larger 
sample of local SFGs from the SDSS and with 
a smaller sample of super compact UV -luminous 
starburst galaxies (jHeckman et al.f2005h placed at 
the same redshift range of the GPs. In particular, 
the position of GPs in the fundamental relation 
between mass and metallicity is discussed in the 
context of their starburst activity and feedback 
processes. 

2. THE DATA 

The sample of galaxies was taken from the list 
of C09. They compiled and analyze d by means 
of SDSS Data Release 7 (SDSS DR7; lYork et al 



20001) imaging and spectroscopy, Hubble Space 
Telescope (HST) optical and GALEX UV imag- 
ing, a well-defined sample of GPs in the redshift 
range 0.11 < z < 0.36 based on color selection cri- 
teria. From their list we selected a sub-sample of 
79 GPs spectroscopically classified as star-forming 
systems (their Table 4), ruling out active galac- 
tic nuclei (AGNs) according to the diagnostic di 



agram [Nil] /Ha versus [Om]/H/3 ([Baldwin et al 
1981). We refer the reader to C09 for a complete 



description of the sample selection. 

We obtained the full set of calibrated one- 
dimensional spectra, i.e., spSpec fits files, from the 
SDSS archive using the SDSS DR7 query toolfl 
The SDSS spectra cover 3800A< A <9200A and 
the fiber has a 3" diameter which, due to the com- 
pactness of the galaxies, include >70% of the GP 
r— band emission^. Emission-line integrated fluxes 
were consistently measured for each galaxy spec- 
trum using the IRAF task splot, by summing all 
the pixels of the emission line after a linear sub- 
traction of the continuum. The flux uncertainties 
are then calculated as 
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N - 



EW 



(1) 



( Perez-Montero fc Diaz 20031) , where a\ is the line 
flux uncertainty, a c is the uncertainty in the posi- 
tion of the continuum, N is the number of pixels 
in the flux measurement, EW is the equivalent 
width of the line and A is the wavelength dis- 
persion. In most cases the measurements are in 
agreement with those based on Gaussian fitting. 

3. RESULTS 

3.1. Green Peas: A Genuine Population of 
Metal-poor Galaxies 

Physical properties and oxygen and nitrogen 
ionic abundances were calculated from the mea- 
sured emission lines using the direct method, 
which is expected to provide more accurate 
abun dances compared to str o ng-lin e methods 
(see Perez-Montero fe Contini 20091 hereafter 



J http: //cas . sdss . org/dr7/en/tools/search/sql . asp and 



http : //das . sdss . org| 
2 We used f fib = io-°- 4 ( m '-,fiber-»v,p.tr.>) ) which is the frac- 
tion of the fiber total light (line+continuum) within the 
fiber aperture 
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PMC09). First, we have estimated the redden- 
ing constant fro m the Ha/H/3 ratio u sing the ex- 



tinction curve of Cardelli et al. | (jl989h . Then, for 



those objects with a reliable measurement of the 
[0 in] A 4363 emission line (~ 70% of the sample), 
we have derived the [O in] electron temperature. 
The [On] electron temperature has been derived 
from it using the relations proposed by PMC09. 
+ , 2+ , and N + ionic abundances have been cal- 
culated using these temperatures and the relative 
intensities of the corresponding brightest emission 
lines. Finally, total O/H and N/O abundance ra- 
tios have been obtained with the usual expressions 
(see PMC09). 

For the full sample, including those galaxies 
with no direct determination of the tempera- 
ture, oxygen total abundances have also been 
estimated from the strong-line calibrator N2 = 
log([Nn] A6584/Ha). We have used the empirical 
relation derived by PMC09, which gives oxygen 
abundances fully consistent with those derived 
from r o -sensitive methods over the whole observed 
range of metallicity. According to these authors, 
N2 has a non-negligible dependence on the N/O 
ratio. Nevertheless, our calculations indicate that 
correction for this dependence in this sample does 
not reduce the dispersion, so we have preferred 
to keep the calculated abundances in the non- 
corrected form. Regarding N/O, this can be de- 
rived using the N2S2 parameter (i.e., [N n]/[S n]), 
minimizing its dependence on reddening or flux 
calibration. 

We have compared the derived properties of the 
GPs with two different samples of SFGs. First, we 
have taken the large sample of emission- line galax- 
ies listed in the Max Planck Institute for Astro- 
physics/Johns Hopkings University (MPA/JHU) 
Data catalog of the SDSS DR 'fj which covers a 
redshift range 0.03 < z < 0.37. We have kept the 
non-duplicated galaxies with a signal-to-noise ra- 
tio (S/N) of at least 5 in all the involved lines. 
Then, we have discarded those objects classified 
as AGNs, using the diagnostic diagram [Nil] /Ha 
versus [Om]/H/3. Finally, we have calculated oxy- 
gen and nitrogen abundances in the same way as 
described above for the sample of GPs. Moreover, 
using the SDSS SFGs with direct estimation of 
the ionic abundances, we have obtained the em- 



03 



30 
25 

20 



^ 15 

1 

5 




This paper C09 




20 



1 b 



10 




7.5 8.0 8.5 9.0 
1 2 + loc(0/H) 



-1.5 -1.0 -0.5 0.0 

log(N/0) 



Fig. 1. — Oxygen and nitrogen abundances. On 
the left, the lined and the no n- lined histograms 
indicate the number of GPs with metallicities cal- 
culated from the direct method and using the N2 
index, respectively. For comparison, metallici- 
ties derived by C09 are shown in the filled his- 
togram. On the right, the lined and the non-lined 
histograms show the corresponding N/O distribu- 
tions. (A color version of this figure is available in 
the online journal.) 

pirical relation between N/O and the N2S2 pa- 
rameter (PMC09) that was consistently applied 
to estimate N/O for all the galaxies used in our 
analysis. Our polynomial fit: 



log(N/0) 



-0.76 



1.94N2S2 + 0.55 (N2S2) 2 
(2) 

where N2S2 = log([Nn]/[S n]), gives a slightly 
higher slope and a much lower dispersion (~0.1 
dex) than the calibration of PMC09. 

In addition, we have taken for comparison a 
sample of 31 nearby (z < 0. 3) UV-luminous star- 
burst galax ies studied by (|Overzier et all 120081 
120091 l2010h . They showed that these galaxies ex- 
hibit strong similarities in their properties with 
the Lyman-break galaxies at high-z, and there- 
fore were proposed as good local Lyman-break 
analogs (LB As). Most massive LB As display a re- 
solved morphology, characterized by the presence 



of a d ominant central object fDCO. IOverzier et al 



3 Available at http://www.mpa-garching.mpg.de/SDSS/ 



2009). In contrast, low-mass LBAs are metal-poor 
galaxies showing close similarities (e.g., morphol- 
ogy, redshift, UV luminosities, size, SSFRs, and 
environment) to GPs (C09). In fact, three LBAs 
are also included in the GP sample. These three 
galaxies are those for which HST optical imaging 
was used to analyze the morphology of the GPs 
(see Figure 7 in C09). To be consistent, the oxygen 
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and nitrogen abundances for the sample of LBAs 
were re-calculated using the line ratios measured 
from SDSS spectra (|Overzier et al.ll2009h . and the 
same empirical calibrations described above. 

In Figure [1] we show histograms with our 
metallicity estimates for the GPs and those de- 
rived by C09 for comparison. The GPs show 
direct metallicities 7.7 < 12+log(0/H)< 8.4, 
with a mean value of 8.05 ± 0.14, i.e., ^20% 
the solar valu^J whereas their measured C(H/3) 
and Ha/H/3 mean values are 0.23 and 3.3, re- 
spectively. Our direct estimates and those de- 
rived using the N2 parameter agree well in ^0.05 
dex. While extinction values are in agreement 
with those obtained by C09, Figure [1] shows a 
clear offset of ^0.65 dex between their metal- 
licity estimates (mean log(0/H) +12 = 8.7) 
and ours. The C09 values were obtained us- 
ing the [NII]A6584 A/[0 II]AA3726, 3729 A ratio 
to estimate log(0/H) + 12 via the theoretical 
calibration based on ph otoionization models by 
Kewlev fc Dopital ( 20021 ). Systematical discrepan- 
cies are usually found in the literature between 
metallicities de rived from theore t ical a nd empir- 
ical methods ( Kewlev fc Ellison 20081 and ref- 
erences therein). We have used the conver sion 
constants provided by Kewlev fc Ellison ( 20081) to 
scale our N2-based oxygen abundance estimates 
to those expected from the theoretical calibration 
(their Table 3). In doing so, we obtained a mean 
difference of 0.23 dex, which is not enough to ex- 
plain the largest offset in O/H shown in Figure [T] 
We conclude that GPs are genuine metal-poor 
galaxies with low extinction values. 

3.2. N/O Ratios and Metallicity 

The nitrogen-to-oxygen ratio is a powerful evo- 
lutionary indicator in galaxies (jPilvugin et al 



20041 : iMolla et al.ll2006l ). We have investigated 



the relation between the N/O ratio and the oxy- 
gen abundance for the sample of GPs, LBAs, and 
the SFGs taken from the SDSS in Figured For the 
SFGs, we found a positive trend with an increas- 
ing scatter toward the metal-rich regime reflecting 
both primary and secondary production of nitro- 
gen in the same metallicity range. At low metallic- 
ities the trend flattens out, possibly a consequence 
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Fig. 2.— N/O vs. O/H. The two-dimensional 
histogram shows the number distribution of SDSS 
SFGs in the logarithmic gray scale. GPs are in- 
dicated by green circles, while LBAs are indi- 
cated with open red squares. LBAs with DCOs 
are marked with open red squares with plus sign, 
whereas those LBAs also belonging to the GP sam- 
ple are indicated with filled red squares. (A color 
version of this figure is available in the online jour- 
nal.) 

of nitrogen being of primary origin, coming mainly 
from massive stars. In the low-to-intermediatc 
metallicity range, most GPs and some LBAs dis- 
play systematically larger N/O ratios for a given 
metallicity compared to the SDSS SFGs. Possi- 
ble reasons for this include extra production of 
primary nitrogen, coming from low-met a llicity 
intermediate-mass s t ars (Renzini fc Volil Il981t 
Gavilan et all l2006t IMolla et all I2006I). pollu- 



tion by Wolf— Rayet stars (e.g., Brinchman n et al 



2008; Lo pez-Sanchez fc Esteban 2010; M onreal-Ibero et al 



20101 ) . or hydrod ynamical effects involving outflow 

and inflow of gas ( van Zee et al .119981: iKoppen" fc Hensler 
20051 ) . This systematically larger N/O ratios could 
explain the difference between our metallicity and 
the higher est imates obtained by C09 using the 
calibration bv lKewlev fc Dopital (|2002l ). The lat- 
ter assumes that N/O is roughly constant be- 
low 12+log(0/H)~8.4, and proportional to O/H 
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above this value. This implies that larger val- 
ues of [Nil] /[On] lead to higher metallicities, 
but this calibration does not take into account 
the dependence of [Nil] /[On] on the variation 
of the N/O ratio at a given O/H. According to 
PMC09, large values of the N/O ratio can enhance 
the [Nil] /[Oil] ratio even in the low metallicity 
regime. 

3.3. Relation Between Stellar Mass, Metal- 
licity and the N/O Ratio 

In FigureEl we studied th e stellar mass— metalli city 
relation (M* - Z, MZR; iLequeux etUI Il979h 
and the relation between N/O and the stellar 
mass (M*-N/0; PMC09) for the GPs, LBAs, 
and SDSS SFGs. Estimates of M* for the three 
galaxy samples have been taken from C09 (GPs), 
Overzier et all (|2009l LBAs), and the MPA-JHU 
catalog (SFGs). They were derived using spectral 
energy distribution fitting and SDSS photometry. 
Thus, stellar masses are expected to to be con- 
sisten t within a typical dispersion of ~ 20 de x 
(e.g., iDrorv et al.ll2004t Ivan der Wei et alJl2006t l. 
This value, is lower than the expected uncertain- 
ties of 0.3 dex quoted for GPs and LBAs as due to 
uncertainties in thei r star formation history (C09; 
Overzier et"alll2009h . 

In the M*-Z plane, SDSS SFGs show a clear 
positive trend which flattens toward higher masses 
and show a relatively large scatter at lower masses, 
in agreement with previous findin gs for local 
galaxies (e.g., iTremonti et al.l l2004i ). Similarly, 
GPs and LBAs increase their metallicity with in- 
creasing masses. Nevertheless, GPs lie more than 
a factor of 2 (> 0.3 dex) below the MZR of SDSS 
SFGs, i.e. at a fixed stellar mass GPs are system- 
atically more metal-poor than normal SFGs. A 
similar shift in the MZR has been observed for the 
less massive Qog(M* ) 



< 



10) LBAs (|Hoopes et al 



20071: lOverzier et alj|2O10l this work) . In terms of 
stellar mass, the observed offset (as large as 1 dex) 
largely exceeds the typical uncertainty quoted for 
the GP masses. 

On the other hand, the SDSS SFGs display 
a correlation between N/O and M*, with higher 
N/O ratios at higher stellar masses, in agreement 
with PMC09 for a different sample of local SFGs. 
The existence of the M+— N/O relation reflects the 
fact that the most massive galaxies evolve more 
quickly and, hence, they should have on average 



higher metallicities and N/O ratios. Though the 
scatter in the M*— N/O relation is large, most GPs 
and low-mass LBAs are roughly consistent with 
the trend of SDSS SFGs, and no systematic offset 
is observed. 

4. DISCUSSION 

The shape of the MZR in galaxies has been 
found to depend on several key processes in galaxy 
evolution, in cluding the efficiency of star for- 
mation (e.g., Leaueux et al. 1979 ;_ Molla fc Diaz! 



20051 : iBrooks et al.ll2007t ICalura et al.ll2009h . the 
action of metal-rich s elective outflows and metal- 
poor gas inflows (e.g ., Larson 1974 : Garnetd 2002: 
Tremonti et~al]|2004l: iFinlator fc Davell2008l) . and 
the p ossible variations in the initial mass func- 
tion (jKoppen et al.l 120071) . Moreover, different 



amounts of dark matter can also assist in some of 



these mechanisms (jDekel fc Silld 119861 ) . 

The GPs follow a relation between mass and 
metallicity that parallels the MZR defined by 
the SDSS SFGs, but is offset >0.3 dex to lower 
metallicities. Interestingly, we find some remark- 
able similarities between the GPs and the pop- 
ulation of [O in] -selected galaxies of similar lu- 
mi nosity in the z rang e 0.29-0.42 recently found 
by ISalzer et al. (2009). These galaxies follow a 
luminosity— metallicity relation that parallels the 
one defined by SFGs, but is offset by a factor 
of more than 10 to lower abundances. On the 
other hand, since GPs and low-mass LBAs lie in 
a similar offset position in the MZR relative to 
normal SFGs, their mentioned similarities appear 
even greater. One possibility to explain their offset 
position in the MZR is that these galaxies could 
be still converting a large amount of their cold gas 
reservoirs into stars. In that case, their low abun- 
dances could be due to their relatively young ages 
compared to normal SFGs. In the same range of 
masses, both GPs and LBAs have much higher 
SSFR (typically >10~ 9 yr _1 ) compared to other 
SFGs of similar mass (C09). Recent studies show 
that galaxies with higher SSFRs or larger half- 
light radii for their stellar mass have systemati- 
cally lower metallic ities (e.g., Tremonti et al.ll2004l : 
Ellison et al.ll2008l ). However, we have found even 
greater under- abundances in the GPs, which have 
high SSFRs but are extremely compact. 

Some models show that in highly concentrated 
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Fig. 3. — M* vs. O/H (left) and M* vs. N/O (rigth). Symbols stand for the same type of objects as 
described in Figure [2] For better comparison the trends observed in the M* — Z plane, green and black lines 
indicating second-order polynomial fits to the GPs and SFGs, respectively, have been overplotted. (A color 
version of this figure is available in the online journal.) 



(typical sizes <3 kpc) low-mass galaxies such as 
the GPs, galactic winds induced by their large 
SSFR are strong enough to escape from their weak 
potential wells, diminishing the observ ed global 
abundances (e.g.. IFinlator fe Dave 2008T). I n con- 
trast, analytical models by DalcantonT i 2007 ) show 
that any subsequent star formation to the out- 
flow will remove their effects on metallicity, un- 
less galaxies have an inefficient star formation. 
Smoothed particle hydrodynamics (SPH) plus TV- 
body simulations have shown that low star for- 
mation efficiencies, regulated by supernova (SN) 
feedback, could be primarily responsible for the 
lower metallicities of low- mass galaxies and their 
overall tr end in the MZR ferooks et al.ll2007h . As 



shown bv lErb et a l. (2006) for SFGs at z ~ 2, the 



constancy in the offset of the MZR suggests the 
presence of selective metal-rich gas loss driven by 
SN winds. 

Inflow of metal-poor gas, either from the out- 
skirts of the galaxy or beyond, can dilute metals 
in the galaxy centers, explaini ng an offset to lower 
abundances in both the MZ R (Mihos & Hcrnquist 



1996[ iBarnes fe Hernquisd H996t IFinlator fc Dave 
20081) and the N/O - O/H diagram. In starburst 
galaxies, a recent cold gas accretion can be due 
to interactions, which eventually increases the gas 
surface density and con sequently the star fo rma- 
tion. As explained by lEllison et al. (2008), the 
dilution of metals due to an inflow can be restored 
by the effects of star formation depending on the 
dilution-to-dynamical timescale ratio. Since this 



ratio depends inversely with galaxy radius, galax- 
ies with smaller radius, such as the GPs, may 
be expected to take longer time to enhance their 
oxygen abundances to the values expected from 
the MZR. In this line, the position of GPs in the 
A/*— N/O relation and the offset observed in the 
N/ O - O/H plane may favor this scenario. Models 
by lK5ppen fc Hensler ( 2005 ) have shown that the 
rapid decrease of the oxygen abundance during an 
episode of massive and rapid accretion of metal- 
poor gas is followed by a slower evolution which 
leads to the closed-box relation, thus forming a 
loop in the N/O - O/H diagram. 

The inflow hypothesis is also strongly sug- 
gested by the disturbed morphologies and close 
companions observed in spatially resolved HST 
images for three G Ps and most LBAs (C09; 



Overzier et al.l [2008J). Recent results revealed 



that galaxies involved in galaxy interactions fall 
>0.2 dex below the MZR of normal galaxies due 
to tidally induced large-scale gas inflow to the 
galaxies' central regi ons (e g iKewlev et al ll2006l 



Michel-Dansac et all 120081 IPeeples et aiT l2009') 



Several iV-body/SPH simulations have shown that 
major interactions drive starbursts and gas inflow 
from the outskirts of the H I progenitor disks (e.g., 
Mihos fc Hernquisd[l99l iRupke et al.ll2010h . also 



supporting this scenario. 

We conclude arguing that recent interaction- 
induced inflow of gas, possibly coupled with a se- 
lective metal-rich gas loss driven by SN winds, may 
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explain our findings and the known galaxy proper- 
ties. Nevertheless, further work is needed to con- 
strain this possible scenario. In particular, future 
assessment of the H I gas properties and the star 
formation efficiency of GPs, as well as the behavior 
of effective yields with mass compared with mod- 
els of chemical evolution, will shed new light on 
the relative importance of the above processes. 

Our results allow us to further constrain the 
evolutionary status of the GPs. These galaxies, as 
well as the low-mass LBAs an d the [O in]-selected 
galaxies bv lSalzer et al.l (J2009J) should be analyzed 
and compared in more detail to elucidate whether 
these rare objects are sharing similar evolutionary 
pathways. Even if this is not the case, their prop- 
erties suggest that these galaxies are snapshots 
of an extreme and short phase in their evolution. 
They therefore offer the opportunity of studying 
in great detail the physical processes involved in 
the triggering and evolution of massive star forma- 
tion and the chemical enrichment processes, under 
physical conditions approaching those in galaxies 
at higher redshifts. 

Forthcoming analysis, based on high S/N 
intermediate/high- resolution spectroscopy and 
deep NIR imaging with the Gran Telescopio Ca- 
narias (GTC), will be used to better constrain the 
evolutionary status of the GPs. 
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